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' Abstract. I have recently combined HST/STIS spectrophotometry with ex- 

\ isting photometric data to analyze the calibration of three standard optical pho- 

. tometry systems: Tycho-2 B^Vt, Stromgren uvhy, and Johnson UBV. In this 

pLn ' contribution I summarize those results, present new ones for 2MASS JHKs, 

and combine them with recent literature results to generate a uniform set of 
zero points for six photometric systems, the above mentioned plus Cousins RI 
and SDSS ugriz. With the exception of the latter system, the zero points use 
, the new Vega spectrum presented at this meeting by Ralph Bohlin. I also discuss 

^ ■ the implementation of these results in CHORIZOS, a Bayesian photometric code 

' that compares multi-filter observational data with spectral energy distributions 

fT^ , to solve the inverse problem of finding the models which are compatible with 

Tij" ' the observations. 

a^ 
o 

o 



1. Introduction 



O \ The inverse photometric problem in astronomy can be defined in the following 

way: one observes a series of magnitudes niohs.p in a series of filter passbands 
, (denoted here by the p index), each one of them defined by a total-system 

dimensionless sensitivity function -Pp(A), and compares them with a family of 
spectral energy distributions (SED, denoted here by the s index) /a,s(A) to find 
^ ' out which of the SEDs is compatible with the observed magnitudes. The first 

^ . required step in that process is to compute the synthetic magnitudes m^^p in the 

r system for each SED in each passband p. For a photon-counting detector, the 
formula is: 

m,p[/,,(A)] = -2.51og,o l^;p^(A)/,_,(A)AdA J + ^^-^ 

A magnitude system r is defined by a reference SED /A,r(A) and a series 
of (relative) zero points ZP^^p for each filter. In principle, one can adjust the 
/A,r(A) in such a way that ZP^^p = 0.0 for all filters, which generates the default 
system for that SED. Examples are given in Table 1. Some magnitude systems 
(e.g. Stromgren-I- Vega) , however, use non-zero ZP^p by definition. 

It is commonly found when one attempts to use a magnitude system with 
zero ZPr,p that, after the data has been fully processed, there are small offsets 
between observed and synthetic values (typically in the hundredths or thou- 
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sands of magnitudes) that have to be corrected^. There are two approaches to 
deal with this issue: the first one, used e.g. for HST photometry, is to leave 
the synthetic magnitudes unchanged and modify mobs,p by changing the calibra- 
tion parameters that are used to convert from instrumental counts to observed 
magnitudes. This leaves the zero points (as defined here) as exactly zero, e.g. 
ZPHST+Vega,p = 0-0 for all filters. Two problems with this approach are that 
the observed magnitudes of a given constant object (e.g. a calibration star) 
change when the reference SED is modified and that it would be hard to apply 
to systems with long histories (e.g. Johnson UBV). An alternative approach, 
which will be followed here, is to calculate the non-zero values of ZP^-.p, leaving 
the published observed magnitudes fixed and introducing changes in the syn- 
thetic ones. This is more practical way to deal with historical systems but it 
means that, for example, the synthetic magnitudes in the Johnson+Vega sys- 
tem^, ?7ijohnson+Vcga,p) where p = U,B, V, are not the same as those calculated 
in the default Vega system (also called VEGAMAG), mvega,p) due to the exis- 
tence of zero points. This one is also a more practical approach for large-scale 
surveys (Tycho-2, SDSS, 2MASS. . . ) where the magnitudes need to be pub- 
lished as soon as possible and it is not until a few years later that accurate zero 
points are calculated. As an example, synthetic SDSS magnitudes, msDSS-i-AB,p, 
where p = u,g, r, i, z, are reported in an AB-bascd system and require a non-zero 
ZPgDSS-i-AB,p in order to be accurately compared with observed magnitudes. 

Table 1. Default magnitude systems {ZPr,p = 0.0 for all p). 

System Reference spectral energy distribution 

Default ST /a,st = 3.63079 • 10"^ erg s"^ cm"^ A"^ 

Default AB /,,'ab = 3.63079 • lO^^o erg s^^ cm^^ Hz"! 

VEGAMAG J x.Vvgii = Vega si")(X'(ruiii 

The relationships between different magnitude systems are easily derived 
from Eqn 1. We start with the magnitudes of Vega in the default ST and AB 
systems, which are given by: 

^ST.(Vega) = -2.51og,o J p^^x)fx,s.XdX ) ' 

, 0.1 /^ /Pp(A)/A,Vega(A)AdA ^ 

n.AB.(Vega) = -2.5 log, [ j p^^,^^,f^^^^/,^,, ) ■ (3) 

"T^ST,p(Vega) and mAB,p(Vega) can then be used to transform from an 
X-|-Vega magnitude system (where X = Johnson, Cousins. . . ) into ST or AB 
magnitudes by applying: 

msT,p = (mx+Vega,p - ZPx+Vega,p) + "^ST,p(Vega), (4) 



^ Given the unintended nature of this effect, some authors refer to these quantities as zero-point 

offsets instead of zero points. 

^That is, the values that should be compared with the observed ones. 
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mAB,p = {rnx+Yega,p - ZPx+Vega,p) + mAB,p(Vega). (5) 

Similarly, if W6 assuniG ZPx-i-Vega,p — Z-Px+ST,;p — Z-Px+AB,p; Can trans- 
form from an X+ST or X+AB magnitude system into an X+Vega system by 
using: 

?^^x+Vega,p = mx+ST,p - msT,p(Vega), (6) 

w^x+Vega,p = 'rnx+AB,p - mAB,p(Vega). (7) 

Some authors list the zero-magnitude f\ or fi, for a given magnitude system 
and filter. Those quantities are easily derived from Eqns. 4 and 5: 

/a,p(0 mag) = /a,ST • 10-°-^[-ST,p(Vega)-ZPx+vega,p] ^ (8) 
/^,p(0 mag) = /^,AB • 10-°-^['"AB,p{Vega)-ZPx+vega,p] . (9) 

It has not been until recently that the quality and stability of the observa- 
tional data and calibrations in astronomy has been good enough to allow for a 
measurement of the photometric zero points of a magnitude system with an ac- 
curacy of < 1%. Such an accuracy is necessary to solve the inverse photometric 
problem without introducing significant systematic errors. In this contribution I 
summarize the recent work that has attained that goal and I generate a uniform 
list of NIR/optical photometric zero points. 

2. A new Vega spectrum 

Vega is the reference BED for many magnitude systems. However, its spectrum 
has been traditionally difficult to measure with high accuracy over large wave- 
length ranges. For example, two of the frequently used SEDs for Vega, those of 
Cohen et al. (1992) [used to link optical and IR] and Bohlin & Gilhland (2004b) 
[optically-based and extended into the IR using a 9550 K Kurucz model\ differ by 
2-3% in flux in the NIR. Part of the difficulty in obtaining an accurate spectrum 
is due to the fact that Vega is a near-pole-on fast rotator (Peterson ct al. 2006; 
Aufdenberg et al. 2006), with a temperature range between 7900 K and 10150 
K over its surface. 

Bohlin (2007) presents in these proceedings a new FUV-to-IR HST/STIS- 

based Vega spectrum that uses a revised CTE correction (Goudfrooij et al. 
2006), which is necessary due to the radiation damage that CCDs suffer in 
space. The new SED uses a 9400 K Kurucz model for long wavelengths that 
improves the overall agreement with the Cohen ct al. (1992) model in the NIR 
to better than 1%. There are still some discrepancies between the two, some 
due to the different absolute calibration for Vega (Megessier 1995 vs. Hayes 
1985) and the rest to the different spectral resolution and version (i.e. year) 
of the used Kurucz models. Nevertheless, the agreement is remarkable (and 
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unprecedented) enough to prompt the usage of this new Vega SED^ as a basis 
for the calculation of accurate zero points. For that reason, I have adopted it 
for CHORIZOS (Mai'z Apellaniz 2004) starting with version 2.1.3 as well as for 
the zero points reported in this contribution. 

3. Results from STIS spectrophotometry 

In two recent papers, Mai'z Apellaniz (2005) and Mai'z Apellaniz (2006), I have 
combined HST/STIS spectrophotometry with existing space-based (for Tycho-2 
BtVt) and ground-based (for Stromgren uvby and Johnson UBV) photometry 
to test the sensitivity curves of those three systems and derive the corresponding 
zero points. Those papers use two samples, the Next Generation Spectral Li- 
brary (NGSL; Gregg ct al. 2004) and the Bohlin (or CALSPEC) sample (Bohlin 
et al. 2001; Bohlin & Gilliland 2004a), which include a large variety of spectral 
types (from O to M), gravities, and metallicities observed with high accuracy 
from 1700 A to 10 200 A. The Vega-based zero points in those papers used the 
Bohlin &; Gilliland (2004b) SED. Here I present an updated summary which 
includes the corrections required to adapt the results to the Bohlin (2007) SED. 

For Tycho-2, observed magnitudes for both Vr and Bt are given in the 
literature, thus allowing a direct test of the validity of the sensitivity curves 
by comparing the observed magnitudes with the synthetic ones derived from 
the STIS spectrophotometry. Also, the Bohlin sample can be used for a direct 
calculation of the magnitude zero points, yielding ZPy^ = 0.036 ± 0.008 and 
ZPbj, = 0.075 ± 0.009. However, as reported in Mai'z Apellaniz (2006), an 
independent measurement for ZPy^, can be obtained from ZPy = 0.026 it 0.008 
(the Bohlin & Gilliland 2004b value, which remains unchanged) and ZPvt-v = 
-0.001 ± 0.005, obtained from the NGSL-I-Bohlin samples, to yield ZPvt = 
0.025 lb 0.010. Combining the two using inverse variance weighting we arrive at 
the final value ZPy^ = 0.032 it 0.006. We can also combine that result with the 
value for ZP Bt-Vt = 0.033 it 0.005 from the NGSL-I- Bohlin samples to obtain 
an independent measurement for ZP^^, = 0.065 ± 0.008. With the two results 
for Bt we arrive at a final value ZPbj- = 0.069 it 0.006. 

For Stromgren u and Johnson U Mai'z Apellaniz (2006), discovered that the 
photometry in the literature was incompatible with the published sensitivity 
curves, thus prompting the derivation of new ones. Given the existing confusion 
between "old-style" (energy-integrating) and "new-style" (photon-counting) sen- 
sitivity curves, Mai'z Apellaniz (2006) also gives the rest of the Stromgren and 
Johnson curves in photon-counting form, thus avoiding the possible errors at the 
hundredth of a magnitude level associated with using the "old-style" functions 
in "new-style" software. It is important to note that neither of those functions 
are "the real sensitivity function" for those filters, especially for the case of 
Johnson U. Instead, what they are is an average over the data published in 
the literature: different observatories (and observers, atmospheric conditions, 
and reduction techniques) are likely to yield different effective sensitivity curves. 



^Available from ftp : //ftp . stsci . edu/ cdbs/ cdbs2/ ciirrent_calspec/alpha_lyr_stis_003 .fits 
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especially to the left of the Balmer jump, where the atmosphere and not the 
telescope/filter/detector system is in most cases the main culprit of the non- 
detection of photons from the source. 

For published Stromgrcn and Johnson photometry, typically magnitudes 
are not given except for V. Instead, colors (b — y, B — V, U — B) or indices {mi, 
ci) arc provided. Therefore, Maiz Apellaniz (2006) only gives color/index zero 
points in order to do a consistent uncertainty calculation (of course, those can be 
transformed a posteriori into magnitude zero points). The random uncertainties 
in the zero points are caused by the finite S/N of the data and the natural 
variation between observing conditions while the systematic errors are caused 
by possible incorrect spectrophotometric or photometric flux calibrations and 
data reductions. The results, adapted to the new Vega spectrum, are given in 
Table 2. 

Table 2. Stromgren and Johnson color/index zero points and associated 
uncertainties/errors using the new Bohlin (2007) Vega SED. 





Stromgren 


Johnson 




b-y 


mi 


Cl 


B-V 


U-B 


zero point 


0.004 


0.157 


1.092 


0.008 


0.021 


random 


0.001 


0.001 


0.002 


0.001 


0.006 


systematic 


0.003 


0.003 


0.004 


0.004 


0.014 



4. Results from NICMOS spectrophotometry 

After Maiz Apellaniz (2006) was published, a number of new spectra were added 
to the CALSPEC sample (Bohlin 2007). In particular there are now 23 stars 
with NICMOS G096-I-G141 grism observations which also have 2MASS JH 
photometry. Of those 23 stars, 14 also have NICMOS G206 observations as 
well as 2MASS Kg photometry. The NICMOS grism observations have been 
corrected of the count-rate dependent non-linearity that affects the detector 
(Bohlin 2007 and references therein). The stars in the sample span a wide range 
of temperatures, from hot white dwarfs to red stars. 

Table 3. 2MASS magnitude zero points using the new Bohhn (2007) Vega SED. 



J 


H 


Ks 


-0.021 ± 0.005 


-F0.009 ± 0.005 


-FO.OOO ± 0.006 



We have applied the technique described in Mai'z Apellaniz (2005) to the 

above sample to measure the zero points for 2MASS JHKg using the new Bohlin 
(2007) Vega SED. Results are shown in Table 3. The comparison between syn- 
thetic and observed magnitudes shows no color terms, indicating that the 2MASS 
sensitivity curves arc well characterized. Also, once the zero points have been 
applied, the distribution of (mobs,p ~ ™r,p)/cp, where Gp is the photometric un- 
certainty and p = J,H,Ks, is closely approximated by a Gaussian of zero mean 
and standard deviation of one. This indicates that the uncertainties are well 
behaved and that the largest contribution to the error budget comes from the 



6 



Ma/z Apellaniz 



photometry, not from the spectrophotometry (as it was the case for the STIS 
analysis in the prevbious section). 

5. Other recent work on photometric zero points 

Other authors have also recently calculated zero points for commonly used op- 
tical/NIR photometric systems. Here I summarize the work of Holberg & Berg- 
eron (2006), Cohen ct al. (2003), and Cohen (2007). They also use a comparison 
of observed photometry with SEDs but with one difference with respect to the 
work in the previous section: their SEDs include not only observed spectropho- 
tometry but also atmospheric models. 

Holberg &: Bergeron (2006) calculate magnitude offsets for Johnson UBV, 
Cousins RI, Stromgren uvby, SDSS ugriz, and 2MASS JHKg. Zero points can 
be immediately calculated by combining them with the observed Vega magni- 
tudes in their Table 1. Their values are derived from a comparison of observed 
photometry of white dwarfs with atmospheric models and use the Bohlin k. 
Gilliland (2004b) Vega STIS + 9400 K Kurucz SED. As we show in table 3 
and in the next section, their results arc in excellent agreement with the filters 
in common with Maiz Apellaniz (2006) with the exception of Stromgren u and 
Johnson U. Those differences are easily explainable by their use of the older 
sensitivity curves for those filters, which Maiz Apellaniz (2006) showed to be 
incorrect (at least in an average sense). 

Cohen et al. (2003) present the calibration of the 2MASS survey and in that 
work they give the zero points for JHKg (which they call zero point offsets or 
ZPOs). It is important to note that their ZPOs are defined with the opposite 
sign as in Eqn. 1 and that they are evaluated with respect to the Cohen et al. 
(1992) Vega SED. In Table 4 we present their results adjusted for those effects. 
Finally, Cohen (2007) presents in these proceedings new ZPs for Tycho-2 BtVt 
derived from their network of FGK dwarfs. Those values are in good agreement 
with the ones in the previous section and are also shown in Table 4. 

6. Uniform zero points 

In order to generate a uniform list of magnitude zero points, I have collected the 
results for the six photometric systems (Johnson UBV, Cousins RI, Stromgren 
uvby, Tycho-2 BtVt, SDSS ugriz, and 2MASS JHKg) mentioned in the previ- 
ous two sections, transformed color /index zero points into magnitude ones where 
required, and adapted them to the new Vega spectrum of Bohlin (2007) where 
needed. All of those systems use Vega as the reference SED with the only excep- 
tion of SDSS. The preferred zero points are given in the sixth column of Table 
4, with alternate values in the eighth column. The third, seventh and ninth 
columns give the references for the sensitivity curve, preferred zero point, and 
alternate zero point, respectively, for each filter. The fourth and fifth columns 
give msT,p(Vcga) and mAB.p(Vcga), respectively. 

Here are some notes on the results in Table 4. The Johnson and Cousins 
values are ultimately linked to the Bohlin Sz Gilliland (2004b) result of ZPy = 
0.026ib0.008. The possible alternate values for Johnson U and Stromgren u from 
Holberg &; Bergeron (2006) are not given because they use a different sensitivity 
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curve, so a direct comparison is not possible. The Cousins RI sensitivity curves 
of Bessell (1983) are in energy-integrating form and should be converted into 
photon-counting form before they are used with e.g. synphot (Laidler et al. 
2005). The preferred Stromgrcn zero points are a combination of the ZPy value 
of Holberg &; Bergeron (2006) with the color /index results of Maiz Apellaniz 
(2006). The alternate SDSS zero points are taken from the SDSS web page and 
they are of lower precision (one hundredth of a magnitude) than the rest. 



Table 4. ST and AB magnitudes for Vega and magnitude zero points 
used by CHORIZOS as of September 2006. 



System 


Filt. 


Ref. 


"iST,p 


W-AB.p 


2Pr,p 


Ref. 




Ref. 








(Vega) 


(Vega) 




alt. 




Johnson-|-Vega 


U 


1 


-0.160 


+0.758 


+0.055 


1,6 








B 


1 


-0.603 


-0.115 


+0.034 


1,6 


+0.027 


7 




V 


1 


+0.005 


+0.003 


+0.026 


6 




7 


Cousins-|-Vega 


R 


2 


+0.547 


+0.184 


+0.030 


7 








I 


2 


+1.208 


+0.419 


+0.017 


7 






Stromgren-|-Vega 


u 


1 


+0.104 


+1.136 


+1.432 


1,7 








V 


1 


-0.765 


-0.140 


+0.179 


1,7 


+0.175 


7 




b 


1 


-0.501 


-0.156 


+0.018 


1,7 


+0.012 


7 




y 


1 


+0.006 


+0.006 


+0.014 


7 






Tycho-2-FVega 


Bx 


3 


-0.650 


-0.072 


+0.069 


1,8 


+0.079 


10 




Vt 


3 


-0.084 


-0.017 


+0.032 


1 


+0.035 


10 


SDSS+AB 


u 


4 


-0.014 


+0.922 


+0.042 


7 


+0.04 


11 




9 


4 


-0.436 


-0.106 


-0.002 


7 


+0.00 


11 




r 


4 


+0.405 


+0.143 


-0.003 


7 


+0.00 


11 




i 


4 


+1.037 


+0.356 


-0.016 


7 


+0.00 


11 




z 


4 


+1.585 


+0.518 


-0.028 


7 


-0.02 


11 


2MASS+Vega 


J 


5 


+2.669 


+0.894 


-0.021 


8 


-0.014 
-0.006 


7 
9 




H 


5 


+3.763 


+1.368 


+0.009 


8 


+0.004 
+0.007 


7 
9 




Ks 


5 


+4.823 


+1.838 


+0.000 


8 


+0.005 
-0.022 


7 
9 



References: [1] Mai'z Apellaniz (2006); [2] Bessell (1983); [3] Bessell (2000); 
[4] http://www.sdss.0rg/dr3/instruments/imager/#filters as of August 2006; 
[5] http : / / www . ipac . caltech . edu/2mass/releases/ allsky/doc/ sec6_4a . html#rsr 
as of August 2006; [6] Bohlin & Gilliland (2004b); [7] Holberg & Berg- 
eron (2006); [8] this work; [9] Cohen et al. (2003); [10] Cohen (2007); [11] 
http://www.sdss.org/dr4/algorithms/fluxcal.html as of August 2006. 

Uncertainties for the zero points are not listed in Table 4 because of the het- 
erogeneous character of the sources: some give them, some do not, and in some 
cases the transformation from color to magnitude uncertainties is not obvious 
given the dependencies between values. In most cases where uncertainties are 
known, they are < 1%. That precision is consistent with the accuracy derived 
from a one-by-one comparison between the preferred and alternate zero points 
in Table 4: all but three differ by 0.010 magnitudes or less and only one by 
more than 0.016 magnitudes. That exception is 2MASS Kg (0.027 magnitudes 
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between the highest and lowest values), which is notorously difficult to calibrate 
using atmosphere models due to the possible existence of very late M- or early 
L-type companions (Holberg & Bergeron 2006). The value derived in this paper 
is located in between the two extremes and is based on observed spectropho- 
tometry, so it should not be affected by such companions, but is derived from 
a small sample (14 objects). For that reason, it would be useful to expand the 
sample of stars observed with NICMOS in order to confirm the result. 

In view of these results, we can conclude that the zero points in Table 4 
allow for a comparison between observed photometry and SED models at a level 
of accuracy of 1% or better. 



7. Implementation in CHORIZOS 

My original motivation when I started working on the issue of zero points was 
to improve the accuracy of CHORIZOS, an IDL code I have worked on for the 
last five years (Mai'z Apellaniz 2004). The purpose of CHORIZOS is to solve the 
inverse photometric problem using multi-filter information with high-precision 
data. Unfortunately, as the attendees of this conference know well, astronomical 
photometry has been often plagued with data whose accuracy does not match its 
precision. In other words, it makes little sense to use photometric uncertainties 
measured in milimagnitudes when the systematic errors introduced by erroneous 
zero points are measured in hundredths of magnitudes. Hopefully, time will 
demonstrate that the zero points presented here do indeed provide the desired 
1% accuracy. They were included in CHORIZOS starting with version 2.1.3 in 
September 2006. Previous versions of the software had more inaccurate values 
(up to 2%), so pasts users may want to consider reprocessing their data. We are 
currently in the process of checking the zero points described in this paper using 
a sample of low-extinction OB stars with multi-filter photometry analyzed with 
CHORIZOS; our preliminary results indicate that they are indeed accurate to 
within 1%. 

In its current (October 2006) version, CHORIZOS calculates synthetic col- 
ors from a grid of model spectral energy distributions and compares them with 
an observed set of magnitudes for an object to select which models are compat- 
ible with the observations. Several stellar and cluster SEDs and 88 filters are 
preinstalled but the user can introduce his/her own. The user may also select 
the type and amount of extinction as two additional parameters and change the 
range and spacing of the SED grid. The multifilter approach has several impor- 
tant advantages with respect to the traditional color-magnitude or color-color 
ones. First, processing all the available information simultaneously allows for 
the easy detection of erroneous data (e.g. an incorrect magnitude) or of objects 
than do not conform to the models (e.g. binary stars in a single-star sample or 
objects outside the assumed parameter range). Also, the combination of multi- 
filter information with a large wavelength coverage makes the simultaneous fit 
of several parameters (temperature, extinction, age, gravity. . . ) possible, thus 
permitting the reduction on the number of a priori assumptions on the proper- 
ties of the object. Finally, since CHORIZOS provides most of its output (text 
and graphical) in parameter space (as opposed to color or magnitude space), it 
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allows for the easy calculation of the properties of the object studied, including 
the possible existence of multiple solutions compatible with the observed data. 

My plans for the future two years include the following modifications to the 
code: 

• A change of the comparison mechanism between synthetic and observed 
data from colors to magnitudes. This will not impact the previously calcu- 
lated results but will allow several other changes, such as adding distance 
or redshift as an extrinsic parameter. 

• Inclusion of spectrophotometric data as input: absolute fluxes, ratios, in- 
dices, and equivalent widths. 

• Inclusion of multicolor indices (Stromgren-like) as input for photometric 
data. 

• Possibility of having detection limits instead of measured magnitudes. 

• Use of arbitrary wavelength gids. 

• Addition of full Bayesian priors. 

• Addition of more graphical capabilities. 

• Inclusion of an algorithm to search and model multiple solutions. 

• Design of a web-based interface and integration with STSDAS synphot. 

In addition, it is possible that small (< 1%) revisions to the zero points will 
be made if new data becomes available. CHORIZOS can be downloaded from 
http : / /www . stsci . edu/~ jmaiz. 

Acknowledgments. I would like to thank Chris Sterken for organizing such 
an opportune meeting and for his understanding during the edition process. I 
also greatly appreciated the useful information exchanges before, during, and 
after the conference with Martin Cohen, Jay Holberg, and Ralph Bohlin (among 
others), which made this contribution possible. 
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Discussion 

Ivezic: All your zcropoint corrections have the same sign. Can you exclude the 
possibility that the STIS calibration is off by ~ 2% in the relevant A range? 

Maiz Apellaniz: It is true that the zeropoint corrections for b — y, B — V,U ~ B, 
and Bt — Vt are all positive (0.004, 0.008, 0.021, and 0.037 magnitudes, respec- 
tively). However, I do not think this is a systematic STIS effect for two reasons: 
First, the uncertainty introduced by the WD calibration (sec talk by Ralph 
Bohlin) is expected to be only 0.002-0.004 magnitudes at most. Second, differ- 
ent authors (M. Cohen, J. Holberg) have recently measured consistent results 
using non-STIS data. 

Sterken: Related to Bessel's remark: the deviations in ci were described by 
us in the LTPV-project as conformity errors. We have used up to 13 different 
uvby systems and we found, empirically, that several such systems were not 
transformable. It would be useful to run your analysis on our results. 

Maiz Apellaniz: I agree that is one possible explanation for the effects I see. 
Indeed, it would be a good idea to test those 13 different systems to check it 
out. 



